Gravitational vacuum polarization as an 
■ alternative to dark matter 

oo ■ 
o 

C\| ■ E. S. Corchero 

■ Departamento de Ffsica. Universidad de Cantabria. Santander. Spain 



O 



February 5, 2008 



O . 

O" 1 Abstract 

^jy It is assumed that the quantum vacuum may be studied as con- 

sisting of two contributions, with positive and negative energy respec- 
tively, which interact but slightly and may be displaced from each 
other. Then it is proposed that dark matter may be just an increase 
CN \ of the quantum vacuum energy, with respect to the normal dark en- 

ergy level, induced by the gravitational field of galaxies or clusters. A 
Q ■ simple model is worked out able to reproduce astronomical observa- 

, tions. 
O . 

An outstanding open problem in present day cosmology is the nature of 
"dark energy" and "dark matter", which toghether contribute more than 
95% of the total mass-energy of the universe. The dark energy corresponds 
^ ■ to a mass density 

Pde ~ HT 26 %/m 3 , (1) 

and the dark matter to about one third this quantity on the averagepQ. A 
common hypothesis is to identify the dark energy with the cosmological con- 
stant introduced by Einstein in 1917. An alternative hypothesis, equivalent 
in practice, is to postulate a vacuum stress-energy tensor, T Miy , of the form 

= pg>", (2) 

g/j-f being the metric tensor, which amounts to assuming an equation of state 
for the vacuum of the form 

p = -p < 0, p = T °, p = -Tl = -Tl = -T|. (3) 
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If the dark energy is due to the quantum vacuum, it is difficult to under- 
stand why its density is not either strictly zero or a density at the Planck 
scale, that is about 10 123 times the observed value, eq.flT]) . We might assume 
that the vacuum contains both positive and negative energy contributions 
which do not cancel completely each other, but then an extremely fine tun- 
ing is requiered in order to get the observed value. In the present paper I 
propose that the dark energy derives from the fact that the positive energy of 
long wavelength components of the vacuum fields is not completely cancelled 
by the other field components whilst its positive pressure is more than com- 
pensated by the negative pressure of other field components. This hypothesis 
reduces substantially the fine tuning needed and, in addition, leads naturally 
to the hypothesis that the vacuum may be "gravitationally polarized" . That 
is, vacuum components contributing positive (negative) energy will tend to 
go to regions with more (less) negative gravitational potential. The purpose 
of the present paper is to show that this assumption may allow interpreting 
"dark matter" as a gravitationally polarized "dark energy" . 

In more detail the argument is as follows. We should assume that all 
short-wavelength components of the quantum vacuum are strongly coupled 
to each other, as shown by the coupling of the corresponding vaccum fluctu- 
ations. For instance, virtual photons with Compton wavelength (~ 10~ 12 m) 
may produce virtual electron-positron pairs and viceversa. This suggests that 
short- wavelengths of the different vacuum fields could not be separated from 
each other and their joint contribution to the vacuum energy, either positive 
or negative, should be not too big in absolute value. In contrast, vacuum 
field components with long wavelength may be almost uncoupled to those 
with short wavelength, but either completely cancel the short-wavelength 
contribution or give a very small, but positive, total vacuum energy. This 
hypothesis offers the possibility that the total vacuum energy and pressure, 
in the absence of gravitational field, fulfils the dark energy eqs.flT]) to (J3]) 
without the need of an extremely fine tuning. 

Now it is natural to believe that the long wavelengths part of the vacuum 
fields may propagate with some independence from the short wavelength 
components. In order to get a simple model, I shall assume that there is 
a cut-of, A c , such that field components with longer wavelength are com- 
pletely decoupled from the remaining vacuum fields. The model so derived is 
clearly too crude and a more plausible model should include some interac- 
tion between the positive energy and the negative energy contributions, thus 
preventing a too big vacuum polarization. This possibility will be considered 
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elesewhere. Here I shall pursue with the development of the simple model by 
fixing the cut-off, A c . It is obvious that it should lie somewhere between the 
Compton and macroscopic wavelengths, but it is difficult to make a precise 
estimate. I shall choose A c ~ lpm, at about the wavelength of visible light, 
which will lead to good agreement with known dark matter properties (see 
below). This value corresponds to a particle mass 

m ~ hc/\ c ~ leV/c 2 



not far from neutrino masses. Thus wavelengths greater than lpm are neg- 
ligible for vacuum fields other than the electromagnetic one. 

The mass density associated to the electromagnetic zero-point field (en- 
ergy l/2hv per normal mode) cut-off at A c is positive and has the value 

Po = S / (nhu)^-du = - I ~4x 10- 1G kg/m 3 . (4) 



cf 
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The essential hypothesis of the model here proposed is that the vacuum 
contribution due to long wavelengths, with density eq.fllj) , and the remaining 
vacuum contributions, with negative density, behave as two non-interacting 
fluids. In particular in a static space-time, the only case to be studied here, 
each part should be in hydrostatic equilibrium in the gravitational field. 

For simplicity I shall consider a static space-time with spherical symmetry, 
where I should study the hydrostatic equilibrium of three non-interacting 
fluids, ruled by 

^ = -&(p J+Pj ),j = 1,2,3, (5) 

where j = 1 corresponds to the long wavelenghts vacuum contribution, j = 2 
to the remaining vacuum fields and j = 3 to the baryonic matter. (Units 
c = G = I will be used throughout this paper.) The gravitational potential, 
is given by 

d$ m + 47rr 3 p f r 2 

<P = — — = — 7, , m(r) = / A-Kplrjr dr, (o) 

dr r 2 — 2mr J 

where p = p x + p 2 + p 3 is the total mass density and p = pi + p 2 + P3 is the 
total pressure at r (all measured in the local frame). It may be realized that 
the sum in j of the three eqs. , combined with eq. flS]) gives the standard 
Tolman-Oppenheimer-Volkoff equation of general relativity, namely 

dp d& , . m + Anr^p . 

/ = -~r {p + p = — 2 5 (p+p)- (7) 

dr dr r A — 2mr 
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It is easy to relate the density and pressure of each component of the 
vacuum with the potential $ if we know the corresponding equation of state. 
For the long-wavelength electromagnetic zero-point field I must assume the 
radiation equation of state, that is 



1 

3' 

which, inserted in the first eq.© , leads to 



1 dpi 
3 dr 



4 

3 H1 



Pi = Pio ex P(- 4$ ) 



(8) 



(9) 



p 10 being the contribution to mass density at infinity (where $ — > 0). As said 
above the contribution of the remaining vacuum fields to both the energy 
density and the pressure should be negative, in order to cancel almost com- 
pletely the positive contribution of the long-wavelength (electromagnetic) 
field. I shall assume that the second contribution, seen as a fluid, is more 
rigid that the first one so that the pressure as a function of the energy den- 
sity is more steep than eq. (jHJ) • As a simple equation of state which allows 
analytical solutions I propose 



(p ) 2 

P2 = 7^^, P20 < °> 
P20 



(10) 



where 7 is a dimensionless parameter and p 20 is the negative contribution to 
the mass density at infinity. Putting this into the first eq. (J5J) I obtain 
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1 d{p 2 ) 
P20 dr 



-<£>' 



P2 + 7 



(ft 



\21 



P20 
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;i + 7 )exp ( 1 - 1 



(11) 

The sum of the two contributions will give the total energy density and 
pressure of the vacuum, namely 



p = p 10 exp(-4$) + 



P20 

7 



1 + 7) exp ( 



;i2) 



V = \pio exp (-4$) + ^0 

3 7 



;i + 7 )exp(--$) -1 



(13) 



Now it is trivial to find values of p 10 , p 2 o an d 7 such that, when $ = 0, we 
get the dark energy values, that is p = —p = p DE with p DE given by eq.(0Q) . 
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I get 
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Pio = (K + l) p DE , p 20 = -Kp DE with K = ^—j » 1. (14) 



Assuming p 10 ~ — p 20 ~ p as given in eq.flU) I obtain K ~ 4 x 10 10 , a very 
high value which means that the density and pressure of the two assumed 
vacuum components cancel quite accurately. 

Eqs.f ll2p and ([13]) may be simplified taking into account that, in the cases 
of interest for us, the potential is very weak, that is |$| << 1. To second 
order in the potential the equations become 



P = Po 



cxp(- !.*) 3 - I cxp ( -i$ 



-2Po $ + yPo $2 + PDE ; (15) 



P = ^Po ex P(~ 4$ ) - 3 Po 



4 1 
-exp(--$) - 1 



Pde - Po $2 - Pde- ( 16 ) 



Typically $ ~ — 10 -6 in galaxies so that, identifying p with the value eq.((4]) I 
obtain — 2p $ ~ 10~ 21 kg /m 3 » Pde- m contrast the dark matter pressure 
predicted by the model is of the same order as the dark energy pressure, 
— Pdei but has opposite sign. The comparison with observations is as follows. 
Recent measurements in the Milky Way[2] give a mass M(< QOkpc) = (4.0 ± 

0.7) x 10 n M© so that the potential is $(60fcpc) M/R = -3.3 x 10~ 7 , 

consistent with the above estimate for $ within the galaxy. On the other 
hand it is known [3 J that the local density, near the solar system, is 0.2 — 
0.5Gev/cm 3 ~ (4 — 9) x 10~ 22 kg/m 3 , in agreement with our model estimate. 
Eq. (Ti~5]) allows writing the "vacuum polarization law" in the simple form 

Vp ~ -2 Po V$ = 2p g, (17) 

g being the gravitational field. But I stress that this relation should have 
limited validity. Indeed, in strong fields the vacuum polarization would suffer 
some saturation so that |Vp| is smaller than predicted by eq. ffTTj) . On the 
other hand for very small potential, $, eqs.f fT5|) and ffl6|) might be modi- 
fied because, being small differences of big quantities, they would be rather 
sensitive to any change in the model. 

We see that the mass density, eq.fTTol) , and pressure, eg-fTTHl) , due to the 
polarized vacuum energy mimic a fluid of cold matter with equation of state 

Pdm = J^Pdm « Pdm, ( 18 ) 
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where p DM = p—p DE is the density above the dark energy level, and similarly 
for the pressure. Thus in the study of dark matter we might use eq. (fTS]) as 
the equation of state of a hypothetical fluid. For instance if eq. (fT8l) is put in 
the hydrostatic equilibrium equation (see eq.©) I get eqs. f|T5|) and (fT6l) to 
leading order in $, that is 

Pdm - -2p $, p DM ^ P $2 - (!9) 

Of course eqs. (|T5|) to (|T8|) are specific for our simple model, in particular they 
rest upon eq. (TTU]) . However the qualitative behaviour of the polarized vacuum 
energy mimicking cold matter is likely valid for any reasonable model, e. g. 
with equation of state different from eq. CTUj) . 

In the following I will study the dark matter problem in a galaxy assumed 
spherical. The study of more realistic galaxies and clusters would be similar 
although more involved. For this study I may neglect the dark energy con- 
tribution, which is very small whenever the dark matter is relevant. In order 
to get the spatial distribution of dark matter I shall solve eq.([6]) . For this 
purpose the Newtonian approximation is good enough, that is 



$' ~ ^ =>• — (rV) = = 47rr 2 (p + p B ) = -8rrr 2 p ® + 4 7 rr 2 / o B , (20) 
r z dr dr 

where p B (r) is the baryonic mass density of the galaxy and I have approxi- 
mated p by the first term of the right hand side of eq. ( 115ft . In the external 
region, where baryonic matter is negligible, the solution of eq. (120]) is 



where 



$= -— sin(pr + <5) ^p = sin (pr + <5) , (21) 
r r 

p = ^/8np G/c 2 , p" 1 = R DM ~ 3 x 10 20 m ~ lOkpc. (22) 

Rdm is the typical radius of the dark matter distribution predicted by the 
model and I have used for p the value (jlj). The solution in the internal region 
cannot be found without knowing the distribution of baryonic mater. As a 
simple model I shall consider a distribution consisting of a constant density, 
p B , within a sphere of radius R , zero outside. Then the regular solution of 
eq.( l20|) for r < R is 

$ = — sin (pr) =^> p = — — sin (pr) — p B . (23) 

2p r r 
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The condition that both $ and $' are continuous at r = R leads to 
Asin(pR) -C sin (/iR + 5) = £*-R=^f^ : 



2p V? 

Acqb(jjlR) -C cos (jiR + 5) = J^- = ^*, (24) 



where eq.( l2"2"j) has been taken into account. These equations allow getting 
two of the parameters A, C and 5 in terms of one of them. Thus I will write 
the dark matter distribution in terms of 5 as follows. Multiplying the first 
eq.(|24"]) times cos (pR) and the second one times sin (pR) and subtracting I 
get 

4tiR 
C sand = — ^-p B 
p 2 



sin (pR) 

- cos (pR) 



^R 3 p B = M B , (25) 



pR 

Mb being the total baryonic mass of the galaxy. In the second eq. (1231) I have 
taken into account that the typical radius of dark matter in galaxies is much 
larger than the radius of baryonic matter , i. e. fiR = Rbaryonic/ Rdm « 1- 
Putting eq. (j25l) into the first eq.(|24l) I obtain 

where I have again neglected terms of order pR. Using the expressions ob- 
tained for A and C in eqs. ( 12T1) and (1231) , respectively, I get the potential, to 
order pR, and the dark matter density in terms of the parameter 6, namely 



3 pR r 2 \ M 



$ ~ — - + 



B 



2 tan5 2R 2 J R ' 
/3 pR r 2 \ M B 

Mr 2M R P n 

$ ~ sin (fir + 5), p ~ ^ sin (pr + 5), r> R. (28) 

r sin o r sin o 

Provided that ir/2 > 5 > pR, as I shall assume (see next paragraph), the 
latter expression shows that, beyond the typical galaxy radius R, the dark 
matter density decreases as 1/r at the beginning and more steeply later 
on, going to zero at r = (n — 5)/p. This behaviour roughly agrees with 
observations. 

It remains 5 as a free parameter which cannot be obtained from our 
model. Thus the normalization (amplitude) of the dark matter distribution 
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is arbitrary and, furthermore, the model predicts bound states of dark matter 
alone, without any amount of real (e. g. baryonic) matter. This possibility 
seems unplausible. Consequently we should assume that such bound states 
are not stable, but stability cannot be studied within our time-independent 
model. Improvements of the model will be studied elsewhere, but in the 
present paper I will simply fix the amplitude of the dark matter distribu- 
tion by comparison with observations. The total dark matter mass may 
be obtained integrating the density eqs.(!27I) and (|28|) in the region where 
p > 0,that is 

Anpr 2 dr / x sin (x + 5) dx ~ M B \ - + 1 ) , 

smdJ^R \sm6 J 

(29) 

where the dark matter mass within the sphere of radius R has been neglected, 
it being of order (pR) 2 (l/tan5 is of order pR, see below). As the dark matter 
mass is typically more than 10 times the baryonic mass in galaxies[3], the 
choice sinS < 0.3 seems appropriate. 

For large r the density eq. (1281) oscillates between positive and negative 
values, which seems unphysical. We cannot use the expedient of substituting 
p = for the predicted density (I28I) when r is larger than the first zero of the 
density function. Indeed this would contradict the relation eq. (1151) between 
dark matter density and potential. Actually when r is large, and therefore 
|$| small, we should take into account several corrections, e. g. deviations 
from sphericity, the influence of neighbour galaxies, or the presence of some 
baryonic matter in the form of gas. I shall study just one correction, namely 
the dark energy term p DE which should be added to p. In order to see whether 
this term makes the density positive everywhere we should take into account 
that the minimum of the density eq. (l2"5|) will happen near pr ~ 3n/2 with 
the value 

2 Po pM B M B 1 

Pmin ~ -J^-JT, cn • x ~ -0.lp pR . 

{Sti/2 — o) smo R smo 

Thus with p given by eq.®, pR ~ 0.1,sin5 ~ 0.3, M B /R ~ 10" 6 , I get 
Pmin ~ —10~ 23 kg/m 3 , far from being cancelled by the dark energy density, 
eq.fll]) . Thus the predicted total density, p + p DE , is still negative in some 
regions, which shows that the model requires some modifications. These 
would be most important for very small densities where the cancellation 
between the p 1 and p 2 is delicate (see eq.ffT2l)). 

A prediction of the model here proposed is the universal value of the 



8 



dark matter radius, estimated in about lOkpc. Thus the extension of 
the region occupied by dark matter is ruled by the parameter //, with the 
result that it will be concentrated in the central region in clusters (whose 
radius is larger than whilst it would extend well beyond the region of 
baryonic matter in galaxies, which roughly agrees with observations [TJ. Of 
course this property is also true if dark matter consists of a gas of particles 
with an equation of state like eq. ([TBI) . 

In summary I have shown that both dark matter and dark energy might 
be interpreted as vacuum mass-energy provided that the vacuum may be 
(slightly) polarized by the presence of gravity. The particular model here 
proposed roughly agrees with observations. It is to be seen whether a more 
sophisticated model may reproduce all the observed properties of dark mat- 
ter. 
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